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SLEEP APNEA SYNDROME DIAGNOSING DEVICE AND SIGNAL 
ANALYZER, AND METHODS THEREOF 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

The present invention relates to a sleep apnea syndrome 
diagnosing device and a signal analyzer, and methods thereof. 

Background Art 

10 It is said that approximately twenty thousand patients with sleep 

apnea syndrome exist in all of Japan. To check whether the patient has 
sleep apnea syndrome, it is essential to attach plural different sensors to 
various parts of a body, such as hands and feet and an abdomen, of the 
patient and perform a test on the patient, as disclosed in Japanese 

15 Patent Laid-open No. Hei 5-200031. 

Further, it is necessary for different specialized doctors to 
analyze test results thus collected, and furthermore it is finally required 
for these different specialized doctors to gather, bringing their analysis 
results together, to draw a conclusion as to whether or not the patient 

20 has sleep apnea syndrome. 

As can be seen from the above, in a related testing method, it is 
necessary to attach plural sensors to a body of a patient, which causes a 
problem that the patient finds preparations for a test burdensome. 
Besides, since the number of necessary sensors is normally 20 or more 

25 and these many sensors need to be attached, the explanation and 
guidance of the test from the doctor to the patient are complicated. 

Moreover, after measurement, plural different specialized doctors 
are required to analyze test results, which causes a problem that a lot of 
time and cost are necessary. In particular, the analysis of the collected 

30 test results requires a lot of time, whereby the development of a method 
capable of temporary screening in a simple manner is desired. Namely, 
the development of a method of temporarily screening many patients as 
to whether they have sleep apnea syndrome, and only when in doubt, 
performing a more detailed test is desired. In particularly, it is thought 

35 that many potential patients exist within Japan, U.S. and so on, whereby 
it is believed that the development of this screening method is urgently 
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needed. 

SUMMARY OF THE INVENTION 
Hence, the present invention is made in view of the 
5 aforementioned problems, and an object of the present invention is to 
provide a sleep apnea syndrome diagnosing device and a signal 
analyzer, and methods thereof. 

In order to accomplish the aforementioned and other objects, 
according to one aspect of the present invention, a sleep apnea 
10 syndrome diagnosing device, comprises: 

a snoring sound collector which collects snoring sound; 
a snoring sound holder which holds the collected snoring sound; 
a correlation coefficient calculator which divides a time axis of the 
snoring sound held in the snoring sound holder into plural cycles and 
15 which sequentially calculates a correlation coefficient between the 
snoring sound of one cycle and the snoring sound of a cycle next to the 
one cycle; and 

an output section which outputs the correlation coefficient 
calculated by the correlation coefficient calculator. 
20 According to another aspect of the present invention, a sleep 

apnea syndrome diagnosing method, comprises the steps of: 

collecting snoring sound and storing the collected snoring sound 
in a snoring sound holder; 

dividing a time axis of the snoring sound held in the snoring 
25 sound holder into plural cycles; 

sequentially calculating a correlation coefficient between the 
snoring sound of one cycle and the snoring sound of a cycle next to the 
one cycle; and 

outputting the calculated correlation coefficient. 
30 According to another aspect of the present invention, a computer 

program product including a medium recording a program for diagnosing 
sleep apnea syndrome, the program being operable to execute the steps 
of: 

collecting snoring sound and storing the collected snoring sound 
35 in a snoring sound holder; 

dividing a time axis of the snoring sound held in the snoring 
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sound holder into plural cycles; 

sequentially calculating a correlation coefficient between the 
snoring sound of one cycle and the snoring sound of a cycle next to the 
one cycle; and 

5 outputting the calculated correlation coefficient. 

According to another aspect of the present invention, a signal 
analyzer, comprises: 

a signal data holder which holds signal data with irregular 
periodicity; 

10 a correlation coefficient calculator which divides a time axis of the 

signal data held in the signal data holder into plural cycles and 
sequentially calculating a correlation coefficient between the signal data 
of one cycle and the signal data of a cycle next to the one cycle; and 

an output section which outputs the correlation coefficient 
15 calculated by the correlation coefficient calculator. 

According to another aspect of the present invention, a signal 
analyzing method, comprises the steps of: 

collecting signal data with irregular periodicity and holding the 
signal data in a signal data holder; 
2 0 dividing a time axis of the signal data held in the signal data 

holder into plural cycles; 

sequentially calculating a correlation coefficient between the 
signal data of one cycle and the signal data of a cycle next to the one 
cycle; and 

25 outputting the calculated correlation coefficient. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a block diagram explaining an example of the 
configuration of a sleep apnea syndrome diagnosing device according to 
30 an embodiment; 

Fig. 2 is a diagram explaining a theory for calculating a 
correlation coefficient between a waveform of snoring sound of one cycle 
and a waveform of snoring sound of the next cycle in this embodiment; 

Fig. 3 is a diagram explaining a method according to this 
35 embodiment for calculating a correlation coefficient on snoring sound 
with irregular periodicity; 
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Fig. 4 is a flowchart explaining an analysis process of snoring 
sound executed by a computer according to this embodiment; 

Fig. 5 is a graph showing a relative sound pressure level of 
snoring sound and its correlation coefficient when regular snoring sound 
5 continues; 

Fig. 6 is a diagram showing a relative sound pressure level of 
snoring sound and its correlation coefficient when hypopnea or apnea 
occurs; and 

Fig. 7 is a diagram showing the configuration of a sleep apnea 
10 syndrome diagnosing device when the analysis process of snoring sound 
according to this embodiment is realized by hardware. 

DETAILED DESCRIPTION OF THE EMBODIMENTS 
Most patients with sleep apnea syndrome snore. Their snoring 
15 sound is regular when they normally breathe, but in hypopnea or apnea, 
the snoring sound is irregular. Therefore, a medical specialist or a 
clinical laboratory technologist catches changes in this snoring sound by 
hearing and judges whether or not a patient has sleep apnea syndrome. 
Hence, this embodiment is intended to analyze the changes in this 
20 snoring sound quantitatively. A more detailed explanation will be given 
below. 

Fig. 1 is a diagram showing an example of the configuration of a 
sleep apnea syndrome diagnosing device according to this embodiment. 
As shown in Fig. 1, the sleep apnea syndrome diagnosing device 
25 according to this embodiment includes a microphone 10, an amplifier 20, 
a low-pass filter 30, an analog to digital converter 40, a digital IO 50, and 
a computer 60. 

The microphone 10 collects snoring sound of a patient who is 
sleeping, and converts the snoring sound to an analog electrical signal. 
30 This analog electrical signal is inputted to the amplifier 20. The analog 
electrical signal of the snoring sound inputted to the amplifier 20 is 
amplified and then inputted to the low-pass filter 30. A high frequency 
band (noise) is eliminated from the analog electrical signal inputted to 
the low-pass filter 30, and then the analog electrical signal is inputted to 
35 the analog to digital converter 40. 

The analog electrical signal of the snoring sound is converted to 
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a digital signal in the analog to digital converter 40, and inputted to the 
digital IO 50. The digital signal inputted from the digital IO 50 is 
inputted to the computer 60, where data on the digital signal is held and 
an analysis of the snoring sound is performed based on the held data. 
5 Then, a result of this analysis is displayed on a display or outputted from 
a printer. 

Next, a snoring sound analyzing method performed by the 
computer 60 will be explained. The snoring sound accompanying 
breathing occurs regularly every time in synchronization with a breathing 

10 cycle. Therefore, after the collected snoring sound is averaged every 
125 ms to find a sound pressure level, a waveform is cut out by a time 
window having a length corresponding to one cycle (about three 
seconds) of breathing of an ordinary person and used as reference data. 
Namely, in this embodiment, 125 ms is a sampling period of the snoring 

15 sound. 

Subsequently, a waveform having one cycle length of breathing 
adjacent to the aforementioned waveform is cut out in the same manner 
and used as comparison data. A correlation coefficient between these 
reference data and comparison data is calculated. Such a computation 

20 is repeatedly performed along a time axis of the snoring sound to 
calculate a result thereof. If the snoring sound is repeated regularly 
accompanying the cycle of breathing, the value of the correlation 
coefficient calculated by such a computation shifts showing a numerical 
value very close to one. If hypopnea or apnea starts and thereby the 

25 snoring sound becomes irregular, the value of the correlation coefficient 
drops sharply at this point in time. Therefore, it is possible to keep 
track of changes in the snoring sound steadily by the value of correlation 
coefficient. 

Note that the snoring sound is a biosignal, and there are 
30 temporal variations in the timing of its occurrence. Hence, even if the 
time axis of the snoring sound is cut out every time by the same time 
window to calculate a correlation coefficient, the value thereof varies. 
To exclude the influence thereof, in this embodiment, a six-second 
reference data moving period and a six-second comparison data moving 
35 period are provided, and correlation coefficients on all combinations are 
calculated while each shifting is performed by 125 ms each time with 



10 



15 



20 



respect to a three-second breathing cycle. Namely, the reference data 
moving period which is longer than one cycle is set, and the comparison 
data moving period which is longer than one cycle is set as well. 

Fig. 3 is a diagram explaining this calculating method in detail. 
As shown in step 1 in Fig. 3, first, a three-second reference period tn to 
ti2 is set from time 0 sec which is a point 0. In this case, the six-second 
reference data moving period is between time 0 sec (Tn) and time 6 sec 



Next, a comparison period t 2 i to t 2 2 is set, being shifted from the 
reference period tn to U2 by 4.5 seconds. In this case, the comparison 
period is between time 4.5 sec (t 2 i) and time 7.5 sec (t 22 ), and the 
comparison data moving period is between time 4.5 sec (T 2 i) and time 
10.5 sec (T 22 ). This comparison period t 2 i to t 22 is set so as to overlap 
with the latter one fourth of the reference data moving period. Namely, 
the reference data moving period and the comparison data moving 
period overlap only by one fourth (1.5 seconds). 

Then, a correlation coefficient between a waveform of the snoring 
sound of the reference period tn to t i2 and a waveform of the snoring 
sound of the comparison period t 2 i to t 22 is calculated. An expression to 
calculate this correlation coefficient is shown by an expression (1). 
Here, it is assumed that there are n pieces of data between time tn and 
time ti 2 , and that there are n pieces of data between time t 2 i and time t 22 . 
In this embodiment, n = 24. Xj shows a data value of the i-th reference 
period, and yr shows a data value of the i-th comparison period, /x 
shows a mean value of data of the reference periods, and /y shows a 
mean value of data of the comparison periods. 



(Ti 2 ). 




... (1) 
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Next, as shown in step 2, in a state where the reference period 
tn to ti 2 is fixed, while the comparison period t 2 i to t 22 is being shifted by 
125 ms each time, respective correlation coefficients between the 
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waveform of the snoring sound of the reference period tn to t i2 and the 
waveform of the snoring sound of the comparison period t 2 i to t 2 2 are 
calculated one by one. The comparison period t 2 i to t 22 is shifted to the 
end of the comparison data moving period. In other words, the 
5 comparison period t 2 i to t 22 is shifted until t 22 = T 22 is attained. 
Accordingly, in this embodiment, 24 correlation coefficients are 
calculated with respect to one reference period tn to ti 2 . Incidentally, in 
this embodiment, the amount of shifting of the comparison period t 2 i to 
t 22 is 125 ms, which corresponds to the sampling period of the snoring 

10 sound, but these two need not necessarily coincide with each other. 

Next, as shown in Fig. 3, the process in step 2 is repeated while 
the reference period tn to ti 2 is being shifted by 125 ms each time. In 
other words, the process in step 2 is performed while the reference 
period tn to ti 2 is being shifted until ti 2 = T i2 is attained. Consequently, 

15 in this embodiment, 24 x 24 = 576 correlation coefficients are finally 
calculated. In step 4, positions of the reference period tn to ti 2 and the 
comparison period t 2 i to t 22 after 576 correlation coefficients are 
calculated are shown. Incidentally, in this embodiment, the amount of 
shifting of the reference period tn to t i2 is 125 ms, which corresponds to 

2 0 the sampling period of the snoring sound, but these two need not 

necessarily coincide with each other. 

Next, a maximum value is extracted from values of the calculated 
576 correlation coefficients and taken as a representative value. This 
representative value is adopted as a value of the correlation coefficient 
25 of the reference period tn to t i2 . Then, both the reference period tn to 
ti 2 and the comparison period t 2 i to t 22 are shifted by three seconds, and 
the process from step 2 to step 4 is repeated. In this embodiment, by 
repeating the aforementioned process, the correlation coefficient 
between the snoring sound of one cycle and the snoring sound of the 

3 0 next cycle on the time axis is calculated sequentially. 

Fig. 4 shows a flowchart when the computer 60 performs the 
aforementioned correlation coefficient calculating process. As shown in 
Fig. 4, first, a variable i is initialized to one, and a variable j is also 
initialized to one (step S10). Subsequently, the three-second reference 
35 period tn to ti 2 is set (step S12), and the three-second comparison 
period t 2 i to t^ is set (step S14). As described above, in this 
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embodiment, the comparison period t 2 i to t 2 2 is set at a position which is 
sifted from the reference period tn to ti2 by 4.5 seconds. 

Next, a correlation coefficient [i, j] between a waveform of the 
snoring sound of the reference period tu to ti 2 and a waveform of the 
5 snoring sound of the comparison period t 2 i to t 22 is calculated (step S16). 
Subsequently, it is judged whether the comparison period t 2 i to t 22 is the 
end of the comparison data moving period (step S18). When the 
comparison period t 2 i to t 22 is not the end of the comparison data moving 
period (step S18: No), the comparison period t 2 i to t 22 is shifted 

10 backward by 125 ms (step S20), one is added to the variable i (step S22), 
and the aforementioned process from step S16 is repeated. 

On the other hand, when the comparison period t 2i to t 22 is the 
end of the comparison data moving period (step S18: Yes), it is judged 
whether the reference period tu to t i2 is the end of the reference data 

15 moving period (step S24). When the reference period tu to t i2 is not 
the end of the reference data moving period (step S24: No), the 
reference period tn to ti 2 is shifted backward by 125 ms (step S26), and 
the comparison period t 2 i to t 22 is returned to the beginning of the 
comparison data moving period (step S28). Then, one is added to the 

2 0 variable j (step S30), the variable i is initialized to one (step S32), and 
the aforementioned process from step S16 is repeated. 

On the other hand, when the reference period tu to ti 2 is the end 
of the reference data moving period (step S24: Yes), a maximum value is 
selected from values of the hitherto calculated correlation coefficients [1, 

25 1] to [i, j] as a representative value of the correlation coefficients (step 
S34). 

Next, the setting of the reference period tn to ti 2 is shifted 
backward from the setting in step S12 by three seconds, and 
simultaneously the setting of the comparison period t 2 i to t 22 is shifted 

30 backward from the setting in step S14 by three seconds (step S36). 
Subsequently, the variable i is initialized to one, simultaneously the 
variable j is initialized to one (step S38), and the aforementioned process 
from step S16 is repeated. 

Fig. 5 and Fig. 6 are graphs each showing an example of a 

35 diagnostic result by the sleep apnea syndrome diagnosing device 
according to this embodiment. Fig. 5 shows a graph of a relative sound 
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pressure level waveform and its correlation coefficient when breathing is 
performed regularly, and Fig. 6 shows a graph of a relative sound 
pressure level waveform and its correlation coefficient when hypopnea or 
apnea occurs. In Fig. 5 and Fig. 6, when the sound pressure level 
5 waveform is calculated, a one octave band-pass filter process (BPF 
process) is performed at a center frequency of 1 kHz with consideration 
given to characteristics of a snoring sound spectrum and an acoustic 
characteristic in a room where the snoring sound is collected. 

As shown in Fig. 5, when the snoring sound is produced regularly 

10 with breathing, the correlation coefficient shifts, showing a value close to 
one. On the other hand, as shown in Fig. 6, when hypopnea or apnea 
occurs, the correlation coefficient shifts, showing an irregular value lower 
than one. Namely, in Fig. 6, as shown by the relative sound pressure 
level waveform, breathing appears to be stopped for 15 seconds from 

15 time 65 sec to time 80 sec, and correlation coefficient values during this 
period are approximately 0.6, which is lower than correlation coefficient 
values at other times. Moreover, it tunes out that correlation coefficient 
values are low at times other than this period, and that breathing is 
irregular. 

20 As described above, according to the sleep apnea syndrome 

diagnosing device of this embodiment, snoring sound of a patient is 
collected by the microphone 10 and analyzed using the computer 60, 
which makes it possible to save time and labor required for a diagnosis 
of sleep apnea syndrome, resulting in cost reduction. More specifically, 

25 regarding the collected snoring sound, a correlation coefficient between 
one cycle and the next cycle is calculated, and thereby it becomes 
possible to express irregularities in the snoring sound numerically and 
perform a diagnosis, which hitherto depends on a medical specialist's 
analysis, by using the computer 60. 

30 Moreover, with respect to the three-second reference period, the 

six-second reference data moving period which is longer than the 
three-second reference period is set, and with respect to the 
three-second comparison period, the six-second comparison data 
moving period which is longer than the three-second comparison period 

35 is set. Simultaneously, correlation coefficients are calculated on all 
combinations of reference periods set by shifting in the reference data 
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moving period by 125 ms each time and comparison periods set by 
shifting in the comparison data moving period by 125 ms each time, and 
a maximum value thereof is adopted as a correlation coefficient of the 
reference period. Accordingly, even in the case of snoring sound with 
5 irregular periodicity, it is possible to appropriately set a time window and 
calculate a correlation coefficient. 

It should be mentioned that the present invention is not limited to 
the aforementioned embodiment, and various changes may be made 
therein. For example, the result of the analysis process by the 

10 computer 60 is outputted in graphical form in the aforementioned 
embodiment, but the output form thereof is not limited to a graph. For 
example, a list of correlation coefficients may be outputted in numerical 
form to a printer, and when the analysis result contains correlation 
coefficients of 0.6 or less, a warning may be displayed on a display of 

15 the computer 60. 

Moreover, the numerical values of time and length used in the 
aforementioned embodiment are all just examples, and the time and 
length are not limited to these numerical values. Further, the reference 
data moving period and the comparison data moving period need not 

2 0 necessarily have the same length, and the amount of shifting of the 
reference period in the reference data moving period and the amount of 
shifting of the comparison period in the comparison data moving period 
need not necessarily have the same length. 

Furthermore, the present invention is applicable not only to a 

25 diagnosis of sleep apnea syndrome but also to diagnoses of diseases 
using other biosounds. For example, by collecting heart sound and 
subjecting data on this heart sound to the aforementioned analysis 
process, diseases such as arrhythmia can be diagnosed. 

Additionally, the application of the present invention is not limited 

30 to signal data acquired from a living human body, and its application to 
an analysis of signal data with irregular periodicity is possible. For 
example, it is also possible to collect data on rotation sound of a gear 
wheel and subject the data to the aforementioned analysis process to 
thereby detect a defect and the like in the gear wheel. 

35 Moreover, the correlation coefficient calculating process is 

realized by software processing of the computer 60 in the 
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aforementioned embodiment, but it can be realized by hardware 
processing. Fig. 7 is a block diagram showing the configuration of a 
sleep apnea diagnosing device when the aforementioned correlation 
coefficient calculating process is realized by hardware processing. 
5 As shown in Fig. 7, snoring sound of a patient is collected by a 

snoring sound collector 100 composed of the microphone 10 and held in 
a snoring sound holder 110. The snoring sound holder 110 is 
composed, for example, of a volatile semiconductor memory device. 

A correlation coefficient calculator 120 divides a time axis of the 

10 snoring sound held in the snoring sound holder 110 into plural cycles and 
calculates a correlation coefficient between the snoring sound of one 
cycle and the snoring sound of a cycle next to the one cycle in sequence. 
An output section 130 outputs the correlation coefficient calculated by 
the correlation coefficient calculator 120. In the aforementioned 

15 embodiment, one cycle is set to about three seconds. 

More specifically, the correlation coefficient calculator 120 
includes a reference data moving period setter 122, a comparison data 
moving period setter 124, a combination calculator 126, and a 
representative value extractor 128. The reference data moving period 

20 setter 122 sets a reference data moving period having a first length 
longer than a length of the aforementioned cycle on the time axis of the 
snoring sound. The comparison data moving period setter 124 sets a 
comparison data moving period having a second length longer than the 
length of the aforementioned cycle on the time axis of the snoring sound, 

25 the comparison data moving period being shifted from the reference data 
moving period by a first predetermined period. In the aforementioned 
embodiment, the reference data moving period having the first length is 
set to six seconds, and the comparison data moving period having the 
second length is also set to six seconds. The first predetermined period 

30 by which the shifting is to be performed is set to 4.5 seconds. 

The combination calculator 126 calculates correlation coefficients 
respectively on combinations of. cycles set by shifting in the reference 
data moving period by a second predetermined period each time and 
cycles set by shifting in the comparison data moving period by a third 

35 predetermined period each time. The representative value extractor 
128 extracts a representative value based on the correlation coefficients 



calculated by the combination calculator 126. 

In this case, the combination calculator 126 calculates correlation 
coefficients on all the combinations of the cycles set by shifting in the 
reference data moving period by the second predetermined period each 
5 time and the cycles set by shifting in the comparison data moving period 
by the third predetermined period each time. 

The representative value extractor 128 extracts a maximum value 
from values of the correlation coefficients calculated by the combination 
calculator 128 as the representative value. 

10 The aforementioned second predetermined period and third 

predetermined period may coincide with a data sampling period of the 
snoring sound held in the snoring sound holder 110. In the 
aforementioned embodiment, both the second predetermined period and 
the third predetermined period are 125 ms. 

15 The output section 130 can also output the correlation 

coefficients calculated by the combination calculator 126 as a graph. 

Besides, as for the analysis process by the computer 60 
explained in the aforementioned embodiment, it is possible to record a 
program to execute this analysis process on a record medium such as a 

2 0 flexible disk, a CD-ROM (Compact Disc-Read Only Memory), a ROM, a 
memory card, or the like and distribute this program in the form of the 
record medium. In this case, the aforementioned embodiment can be 
realized by making the computer 60 read the record medium on which 
this program is recorded and execute this program. 

25 Moreover, the computer 60 sometimes has other programs such 

as an operating system, other application programs, and the like. In 
this case, by using these other programs in the computer 60, a command, 
which calls a program to realize a process equal to that in the 
aforementioned embodiment out of programs in the computer 60, may be 

30 recorded on the record medium. 

Further, such a program can be distributed not in the form of the 
record medium but in the form of a carrier wave via a network. The 
program transmitted in the form of the carrier wave over the network is 
incorporated in the computer 60, and the aforementioned embodiment 

35 can be realized by executing this program. 

Furthermore, when being recorded on the record medium or 



13 



transmitted as the carrier wave over the network, the program is 
sometimes encrypted or compressed. In this case, the computer 60 
which has read the program from the record medium or the carrier wave 
needs to execute the program after decrypting or expanding it. 



